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Abstract 

A  generalized  elgenfunction  expansion  method,  Churchill's 
method,  and  a  transform  method  are  used  to  investigate  the  unique- 
ness of  the  solution  of  the  equation 


■^  p(y)  ^  <t>(x,y)  +q(y)(|)  +r(y)  — ^  =  0  ,     -00  <  X  <  00  , 


hx 


subject  to  the  boundary  conditions 


and 


^^{x,0)   +ao<t>xx(x,0)  +P^(})(x,0)  =  0 


(})y(x,l)  +a^(t)^(x,l)  +p^(j)(x,l)  =  0 


0  <  y  <  1 


11 


1 .  Introduction 

A.  Welnsteln  [1]  showed  that  if  <j)(x,y)  is  a  potential  function 
which  is  required  to  satisfy 


(1-1)  ^xx^^'^^   ^^^yy^^'^^  =  0  ,  -00  <  X  <  CD  , 


0  <  y  <  1  , 


(1.2)  'j)y(x,0)  -  0  , 


(1-3)  <t'y(x,l)  =  P<|)(x,l)  ,  P  >  0  , 


and  if  A   is  the  unique  positive  root  of 


(1.4)  JT*   tanh  /v^  =  p 
where  p  is  a  constant,  then 

(1.5)  Hx,y)  =  [Aq  cos  X  fh^  +   Bq  sin  x  fh^  ]   cosh  y  v/V* 

is  the  only  bounded  function  which  satisfied  the  above  conditions. 

Weinstein's  proof  of  this  is  based  on  a  completeness  theorem 

connected  with  an  eigenvalue  problem  which  he  introduces  in  the 

following  way.   Let  <t)^^(x,y)  in  (1.1)  be  replaced  by  -A^(y)  while 

<|)   ,  ij)  and  ^   in  (l.l)-(1.5)  are  respectively  replaced  by  V'   »  l/' 
y  y        y  yy    y 

and  tp .      These  substitutions  impose  the  conditions 


(1.6)  '^yy^y^  -Miy)   =  0  ,  0  <  y  <  1  , 


(1.7)  ^y{0)    =    0    , 

(1.8)  ^y(l)  =  pti'(l)  ,  P  >  0  , 


and  these  conditions  define  a  standard  Sturm-Llouvllle  eigenvalue 
problem.   The  eigenvalues  are  the  numbers  which  satisfy 

(1.9)  /^  tanh  /A  =  p  ,  p  >  0  . 

The  eigenvalues  are  all  real  and  we  can  consider  them  as  ordered 
with  respect  to  absolute  magnitude.   The  sole  positive  eigenvalue 
will  be  denoted  by  A*.   The  eigenfunctions  are  given  by 

(1.10)  Tp^iy)   =  \x^   cosh  Y  yir  ,  n=  0,1,2...  , 

and  they  satisfy  the  orthogonality  condition 

1 
f    i^^iy)   tni(y)dy  =  0  ,  m  ^  n  . 

0 

It  is  well  known  that  the  set  (l.lO)  is  complete.   This  Implies 
that  the  function  ^{x,y)    can  be  uniquely  expressed  in  the  form 

oo    Vn^y^ 
<t>(x,y)  =  YZ. — 1 • 


n=0   ,■   p 

^^(y)dy 

0 


In  this  expansion  the  coefficient  h  (x)  is 


^^{^)   =  ^  Hx,y)^^( 


y)dy 


and  It  must  satisfy 


^n^^^  =,/   ^xx^^'y^^n^y^'^y  =  -/  <t'yy(^'y)^n^y)^y 


0 


0 


^y^n  ^^ny 


+  /   Hx,y)^yy(x,y)dy 
0   0 


-  A  h  (x) 


That  Is,  we  must  have 


h^{x)  =  a^  cos  X  /X^  +  b^  sin  x  /A^ 


Hence  if  ^{x,y)    is  to  be  bounded  in  the  strip,  0  <  y  *=  1, 

-oo  <  X  <   CD  ,    then  the  only  non  zero  coefficient,  h  ,  that  can  be 

admitted  is 


hj^  =  h*  =  a*  cos  X  Ta*  +  b*  sin  x  J\* 


and  this  leads  to  the  result  (1.5). 

Let  us  turn  now  to  the  equation 


(1.11) 


<t'xx(^'y^+  <l>yy(x,y)  =  0  , 


-00  <  X  <  OD   , 

0  <  y  <  1  , 


subject  to  the  boundary  conditions 


(1.12)  (t5y(x,0)  =  0  , 
and 

(1.13)  <j)y(x,l)  +a(t)^^(x,l)  =  0  ,  a  >  0  , 

where  a  Is  a  positive  constant.   If  we  attempt  to  analyze  the  solu- 
tions of  this  system  by  following  Welnsteln's  elgenf unction  method 
we  are  led  to  the  eigenvalue  problem  defined  by 

(1.14)  ^yy(y)-At(y)  =  0  ,  0  <  y  <  1  , 

(1.15)  ^^(0)  =  0  , 

(1.16)  ty(l)  =  aAt(l)  . 

This  Is  not  a  standard  Sturm-Llouvllle  problem  because  the  elgen- 
parameter  appears  In  the  boundary  condition  (l.l6).   The  Implica- 
tions of  this  are  substantially  different  from  those  of  the  eigen- 
value problem  discussed  above.   In  order  to  see  this,  take  a  =  1. 
For  this  value  the  eigenvalues  must  satisfy 

tanh  /a  =  /a  . 

Except  for  A  =  A  =  0  the  eigenvalues  are  negative  and  we  consider 
them  to  be  ordered  with  respect  to  absolute  magnitude.  The  eigen- 
functions  {^  (y)}  are 

(1.17)  ^^(y)  =  U^  cosh  y  /r   ,  n  =  0,1,2...  ; 


and  they  must  satisfy  the  generalized  orthogonality  condition 


± 
(1.18)       ^^(1)^^(1)  -  f    ^^(y)^j^(y)dy  =  0  , 


m  ^  n 


0 


Now  if  we  wish  to  proceed  in  accordance  with  Weinstein's  method  we 
need  to  know  whether  or  not  an  arbitrary  twice  dif ferentiable 
function  f(y)  can  be  expressed  in  the  form 


(1.19) 


00  oo 

f (y)  =  \         k  ^  (y)  =  ^  +  "S    {,   ii    (y) 
n=0  n=l 


where  the  ^^^ '  s  are  given  by  (1.17).  If  the  representation  (1.19) 
is  valid  with  the  series  uniformly  convergent  for  0  ■<  y  <  1,  then 
the  coefficients  I     are  fixed  by 


^^•20^    ^n  V.g  __.lg 


IJL   sinh  Va 
^n      ^  n  ■- 


X 

^^(l)f(l)  -J  f(y)^^(y)dy 


0 


This  suggests  the  tentative  association  of  f(y)  with  the  series 
determined  by  (1.20),  that  is. 


oo 


(1.21)   f(y)  -v  i  +2  ^ 


i, 


^(i)f(i)  -  ^    f(y)^n(y)dy 


0 


n=T 


^     2  > 

l-L   sinh  Jl\ 
n       ^  n 


^^(y) 


This  association,  however,  must  be  rejected  because  if  we  choose 
f{y)  =  y  and  note  that 


^  (1) 


1)  -  r  y^. 

0 


y  •'^ri^y^^y  =  °  ' 


n  >  1  , 


2 
we  see  that  (1.21)  forces  us  to  associate  y  with  the  constant  H. 

This  shows  that  the  set  of  eigenf unctions  (1.17)  is  inadequate  for 

expansion  purposes.   We  therefore  conclude  from  the  foregoing 

remarks  that  the  applicability  of  Weinstein's  method  to  the  case 

in  hand  depends  on  finding  a  set  of  functions  (x  (y)}  which  contains 


^n 


ii'   (y)}  and  admits  the  expansion 


00 

f(y)  =  IZ  Cn>^n^y^  • 
n=0 


The  major  part  of  this  report  is  concerned  with  the  develop- 
ment of  expansions  which  allow  an  extension  of  Weinstein's  method. 
In  Section  2,    using  a  Polncare-Birkhof f  formulation,  we  show  how 
the  elgenfunction  method  can  be  applied  to  the  equation 


(1.22)   ^  p(y)  A  (t)(x,y)  +  q(y )<{)  +  r(y )  ^  =  0  , 
ay     5y  ^^^ 


-00  -=  X  <  oo  , 

0  <  y  <  1   , 


subject  to  the  boundary  conditions 

(1.23)  (l)y(x,0)  +a^(t)^^(x,0)  +pQ<})(x,0)  =  0  ; 
and 

(1.24)  ())y(x,l)  +a^(j)^^(x,l)  +  p^(t)(x,l)  =0  . 

Various  analyses  of  the  system  ( 1. 22)-( 1. 24)  appear  in  the 
literature  but  as  far  as  the  author  knows  these  analyses  depend  on 
imposing  restrictions  on  the  real  a ' s  and  P's.   One  of  the  features 
of  the  development  below  is  that  these  constants  are  left 
unrestricted. 


Partial  differential  systems  of  the  above  type  arise  of 
course  in  connection  with  uniqueness  theorems  for  the  associated 
non-homogeneous  equations.   In  many  applications,  however,  they 
also  arise  in  just  the  above  homogeneous  form.   For  example,  the 
the  system  (l . 22) -( 1. 24 )  emerges  for  consideration  in  hydrodynamics 
in  connection  with  flows  in  a  horizontal  channel  with  rectangular 
cross  section.   A  basic  problem  is  to  discover  relations  between 
certain  parameters  which  insure  a  parallel  flow.   Problems  of  this 
type  led  to  the  analysis  contained  in  this  report.   An  example  is 
discussed  in  Section  2. 

Sections  3  and  4  are  devoted  to  the  analysis  of  (1. 22)-{ 1 . 24 ) 
by  methods  related  to,  but  different  from  the  method  of  Section  2. 
Section  3  contains  an  extension  of  a  method  devised  by  R.  Churchill 
whereby,  under  certain  circumstances,  an  eigenvalue  problem  with 
boundary  conditions  dependent  on  the  eigenparameter  can  be  reduced 
to  a  standard  problem.   Section  4  contains  a  discussion  of  the  use 
of  a  transform  method  for  the  study  of  uniqueness  questions  about 
(1.22)-{1.24). 


2.  A  Generalized  Eigenfunction  Expansion  and  the  Extension  of 
Welnstein's  Method 

The  equation 

;s  N  >2, 

(2.1)   ^  p(y)  ^  Mx,y)  + qCy)-!) +r(y)  ^=0,      -aD<x<oo, 

0  <  y  <  1   , 
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subject  to  the  boundary  conditions 


(2.2)  ^^{x,0)  +aj^^{x,0)+^jix,0)    =    0  ; 


(2.3)  (i)y(x,l)  +a^<t)^^(x,l)  +p-^(f)(x,l)  =  0  , 

possesses  the  trivial  solution  (})  =  0.   Are  there  any  other  solutions 
which  are  bounded  in  the  strip  -oo  <x<oo;0^y_^l?   We  assume 
that  p,  q  and  r  are  real  functions  of  y  such  that  p(y)  >   0  and 
r(y)  >  0  for  0  j^  y  ^  1.   We  also  assume  that  the  a's  and  P's  are 
real  but  otherwise  unrestricted. 

In  order  to  analyze  (2.l)-(2.3)  let  us  Introduce  the  eigen- 
value problem  defined  by  the  equations 

(2.4)  ^  p(y)  A  ^{y)   +q(y)^  -Ar(y)V/  =   0    ,  0  <  y  <  1  , 


(2.5)  ^y(O)  +pQ^(0)  =  a^-KfiO)    ; 


(2.6)  ^  (1)+  P-L^d)  =  a^^^(l)  . 


This  problem  can  be  regarded  as  one  which  is  suggested  by  the 
application  of  the  method  of  separation  of  variables  to  (2.1)-(2.3). 

The  elgenfunction  ip    (y)  ^  0  is  a  non-trivial  solution  of 
(2.4)-(2.6)  which  corresponds  to  the  eigenvalue  A  =  A  .   Such  a 
function  is  unique  to  within  a  multiplicative  factor.   These 
eigenfunctions  constitute  a  denumerable  set  {^'^^(y))  such  that 
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1 

(2.7)   /   r(y)^^(y)^^(y)dy-p(l)a^^^(l)^^(l)  + p(0)a^^^(0)^^{0)  =  0  , 
0 

m  ^   n    . 


The  denumerable  set  of  eigenvalues  {A  }  has  the  point  at  infinity 
as  its  only  limit  point;  and  we  take  the  eigenvalues  as  ordered  in 
such  a  way  that 

l^o'  -  '^l'  -  1^2'  •••  1   \\\    1   •••  • 

The  sequence  of  operations  defined  by  the  left-hand  side  of 

(2.7)  occurs  so  frequently  in  the  sequel  that  it  is  convenient  to 
use  a  special  symbol  for  it.   We  define  Q[n,x]  by 

1 
J      r(y)0(y)x(y)dy 
'0 

(2.8)  Q[0(y),x(y)]  =    -  p(l)aiO(l)x(l) 

+  p(0)a^O(0)x(0) 

We  are  interested  here  in  the  problem  of  representing  a  twice 
differentiable  function  f(y)  in  the  form 

00 
n=0 

where  the  set  (x  (y))  contains  the  set  {^^(y)}.   A  basic  expansion 
formula  can  be  obtained  from  a  study  of  the  Green's  functions 
associated  with  (2.4)-(2.6).   This  function  G(y,Ti,A)  is  defined  by 
the  equation 


(2.9) 


LG(y,Ti,A)  -  Ar(y)G  =  5(y-  r])    , 


0  <  y,T]  <  1 


subject  to  the  boundary  conditions 


(2.10)  Gy(0,Ti,A)  +p^G(0,Ti,A)  =  cx^AG(0,Ti,A)  ; 
and 

(2.11)  Gy(l,Ti,A)  +p^G(l,T],A)  =  a^AG(l,Ti,A)  . 


The  symbol  L  is  used  to  denote  the  operation  defined  by 


(2.12) 


LO  ^  ^  P(y)  ^  0(y)  +  q(y)0(y)  ; 


and  with  respect  to  this  we  have  the  fundamental  and  indispensable 
identity 

t  t    t 


(2.13) 


Xl^dy 


p(?(Py-Xyn) 


+  /   OLxdy 


The  symbol  5{y-Ti)  on  the  right-hand  side  of  (2.9)  denotes  the 
generalized  function  such  that  when  it  is  multiplied  by  a  piecewise 
continuous  function  F(y)  (with  a  finite  number  of  finite  jumps) 
and  integrated,  the  result  means 

b 

/  5(y- Ti)F(y)dy  =  I  [F(ti  +  0)+  F(ti  -0)] 
a 

when  a  <  T)  <  b. 

The  function  G(y,'n,A)  can  be  expressed  in  the  form 
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where  each  of  0{u,r\,'K)    and  a3(A)  is  an  entire  function  of  A  =  z 
regarded  as  a  complex  variable.   The  zeros  of  cd(A)  are  just  the 
eigenvalues  of  (2.4)-(2.6).   A  complex  variable  method  depending 
on  the  use  of  the  Cauchy  integral  formula  coupled  with  the  theory 
of  residues,  can  be  used  to  show  that  if  C  is  a  circle  of  radius 
p  centered  at  the  origin  of  the  complex  (A=z) -plane  and  containing 
the  first  m+1  eigenvalues  A  ,  n  =  0,l,2,...,m,  then 

(2.14)  lim  ^    '  fi  GiZ^k^  dz 

p-^oo   ^^  '^ 

is  a  null  function  while 


(2.15)  llm  1       u&  M^iIL^  dz  =   0 

p-^oo    27ri    y 

and 

(2.16)  llm       ^     §r   G^1>J]>^)    dz   =   0    . 

p->oo  ^ 

1 
These   results   imply   that   if  F(y)    is    such   that      /      F    (y)dy  exists, 

then  ° 

1 
/     F(y)G(y,Ti,z) 

(2.17)  lim       ^    #   -^ dy  =   0    . 


^    27ri 
p — >  oo 


C 


Also   if   f(y)    is   a   twice   dif f erentiable   function   such   that 


1] 


Lf(y)  =  F(y) 
then 


(2.18)     f(Ti)  =  -  lim   2^  ^  Q[G(y,Ti,z),f(y)]dz  . 

p^  oo 


The  formula  (2,l8)  is  called  the  Poincare-Blrkhoff  formula. 
It  was  noted  in  a  less  general  form  by  Poincare  [2]  during  some 
work  on  a  special  problem  in  partial  differential  equations.   Then 
Birkhoff  [5]  proved  the  formula  for  an  ordinary  n   order  boundary 
value  problem  subject  to  certain  regularity  assumptions  and  with 
the  eigenparameter  absent  from  the  boundary  conditions.   Later, 
Tamarkin  [4]  showed  that  the  formula  is  valid  for  a  wide  class  of 
boundary  value  problems  with  the  parameter  present  in  the  boundary 
conditions.   Since  then,  the  formula  has  been  proved  by  Wilder  [5], 
Langer  [6],  Rasulov  [7]  and  others  under  less  restrictive  conditions 
than  those  used  by  Tamarkin.   The  proofs  of  (2.l8),  as  given  by  the 
authors  noted  above,  depend  upon  explicit  asymptotic  estimates  of 
the  behavior  of  eigenfunctions  and  eigenvalues  as  A  — »•  00  .   For  a 
proof  of  (2.18)  with  respect  to  the  second  order  system  (2,9)- 
(2.11);  and  one  which  does  not  depend  on  specific  asymptotic  evalua- 
tions, see  Peters  [8]. 

The  formula  (2.l8)  leads  to  the  expansion  of  f(y)  into  an 
infinite  sum  of  residues.   If  C   is  a  circle  with  center  at  A^ 
containing  no  other  eigenvalue  we  have 
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(2.19) 


00   -,    r 
f(y)  =  -) p^   r  Q[G(t,y,z),f(t)] 


n=0 


dz 


n 


00 

5 


Q 


^  7^  G(t,y,z)clz,f(t) 


n 


If  (X»(z)  has  a  zero  of  order  k  at  z  =  A  ,  the  corresponding  term  in 
the  expansion  (2.19)  is 


(2.20) 


Q[?^(t,y,n),f(t)] 


where  9^(t,y,n)/(z  ->^„)  comes  from  the  Laurent  expansion  of 
G(t,y,z)  for  the  neighborhood  of  z  =  A  .   The  function  9      can  be 
obtained  by  substituting  the  expansion 


0 


00 


G(t,y,z)  =  y—   (z  -  A^)*  0_^(t,y,n)  +  T~   [z  -   A^)*  9^(t,y,n) 


in  the  equation 


LG(t,y,z)  -A^r(t)G  -  (z  -A^)rG  =  5(t  -  y) 


and  the  boundary  conditions 


G^(0,y,z)  +PQG(0,y,z)  =   a^A^G(0,y , z)  + a^(z  -  A^)G(0,y,z)  , 


G^(l,y,z) +  p^G(l,y,z)  =  a^A^G(l,y,z)  + a^(z  -A^)G(l,y,z) 


which  define  G(t,y,z).   We  find  from  these  equations,  after 
equating  coefficients  of  like  powers  of  (z  -A  ),  that  the  circum- 
flexed  quantities  must  satisfy 


(2.21)  \{t,y,n)  =  7^^^(t)  , 

(2.22)  (L-  A^r)e^  =  re^  +  5(t  -y)  , 

(2.23)  (L-V^'^j  =  ^^J+1  '  ^    =  l,2,...,(k-l)  , 


with  the  boundary  conditions 


(2.24) 

9j^(l,y,n)  +P-^9j(l,y,n)  =  a^A^Q^  (l,y,n)  +  a^9j^^(l,y,n)  ; 


to  be  satisfied  for 

J  =  0,1,2,. ..,(k-l)  . 

a* 

The  above  equations  can  be  satisfied  only  if  the  functions  9.(t,y) 

satisfy  the  compatbility  conditions 

(2.25)  <^[f^{t),'e.^^{t,y,n)]    =   0 

for 

j  =  l,2,...,(k-l)  ; 

and 

(2.26)  -^[f^{t),e^{t,Y,n)]    =  ^^(y)  . 
If  z  =  A  is  a  simple  zero  of  a3(z)  we  find  from  (2.21)  and 


n 


(2.26)  that 
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0-L(t,y,n)  =  y^f^it) 


■'i[i'^it),7j^{t)]    =  ^^(y) 


^  (y) 

7n  =  - 


^       <^[f^{t),f^{t)] 


and  therefore  (2.20)  Is 


Q[^„(t),f{t)]V'^(y) 
Q[0  (t,y,n),f(t)  =  -    "^  Jl 


n  n 


If  each  zero  of  03(2)  is  simple,  the  expansion  (2.19)  becomes 

00  Q[t^(t),f(t)]^^(y) 


(2.27)  f(y)  =  5 


n'  ' '  '  '  -•  ^n' 


feu   Q[^^(t),^^(t)] 

If  z  =  A  Is  a  zero  of  cd(z)  with  multiplicity  k  ^  1  we  can 
see  from  (2.23)  that  9, (t,y,n)  Is  not  an  elgenfunctlon  but  must 
satisfy 

(2.28)         [i  (L-A^r)}^""^  0^(t,y,n)  =  7^f^{t)    . 

For  this  case  9-,(t,y,n)  is  called  a  generalized  elgenfunctlon. 

We  are  now  in  a  position  to  apply  the  foregoing  results  to 
the  problem  stated  at  the  beginning  of  this  section.   For  any  finite 
X  the  function  <})(x,y)  can  be  expanded  in  the  form 


00 
(2.29)        (t>(x,y)  =  -  ^    Q[0  (t,y,n),.{)(x,t)] 

n=0 
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The  boundedness  of  (|)(x,y)  depends  on  the  boundedness  of  each  term 
(2.50)  (7n(x)  =  Q[0^(t,y,n),(t)(x,t)] 


a 


s  a  function  of  x.   If  we  differentiate  (2.30)  with  respect  to  x 


and  express  the  derivative  of  the  right-hand  side  In  terms  of  6 
using  (2.23),  we  find 

1 


^n^""^  =  "/  ^i(t,y,n)L(|)dt 


0 
-  p(l)a^<l)^(x, 1)9^(1, y,n) 

+  p(0)a^(t)^^(x, 0)9-^(0, y,n) 
which  reduces  to 

We  can  see  from  the  last  result  that  If 

-  dx 


then  a    (x)  must  satisfy 


(2.31)  [D^  +  A^]^  cTn(x)  =  0 


where  k  Is  the  order  of  the  pole  of  G(t,y,n)  at  z  =  A  .   The 


general  solution  of  (2.3I)  Is 


16 


lx/T~   k-l 
(2.32)  a^(x)  =.  e    "^  ^  a^  x 


-Ix/A   k-l 
+  e     "^  5    b  .x-^  . 

Therefore  if  rr  (x)  is  to  be  bounded  for  all  x  we  must  choose 


a^j  =  b^.  =  0  ,    j   =   1,2,..., (k-l) 


Then  what  remains,  namely 


ix/A^       -ix/tT^ 

_.-.    „   _      +  b   ( 
n'  '    no  no 


a^{x)   =  a^^e    '^  +  b...e   '  ^ 


Is  bounded  if  A   is  a  real  non-negative  eigenvalue.   If  A  does  not 
n  DO  j^ 

satisfy  this  condition  then  we  must  also  take 


a   =  b   =0 
no    no 

in  order  to  have  a   (x)  bounded. 

We  conclude  that  the  answer  to  the  question  about  the 
existence  of  bounded  solutions  of  (2.1)-(2.3),  other  than  <{)(x)  =  0, 
depends  on  the  disposition  of  the  eigenvalues  of  (2.4)-2.6).   For 
example,  if  the  parameters  of  (2.4)-(2.6)  are  such  that  all  of  the 
eigenvalues  are  negative  then  the  only  bounded  solution  of  (2.1)- 
(2.4)  is  ({)(x)  =  0. 

The  method  explained  in  this  section  can  be  applied  to  the 
two-dimensional  linear  analysis  of  the  flow  of  a  gravitating 
incompressible  and  inviscid  liquid  whi  r-h  is  confined  to  an 
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Infinitely  long  open  channel  with  a  horizontal  bottom  and  vertical 
walls.   Let  the  density  p  of  the  liquid  be  constant;  let  the  only 
body  force  be  the  gravitational  force  pg  acting  In  the  direction  of 
the  negative  y-axls  and  let  the  x-axls  coincide  with  the  bottom  of 
the  channel.   The  basic  hydrodynamlcal  equations  can  be  written  In 
the  dlmenslonless  form 


(2.33) 


Su-    SUo 

^  +  ^  =  0  , 


(2.3^) 


r    -^s  Su,       Su, 

^  |ii  +  u^    1  +  u^  ^ 
dx    -L  dx     ^  dy 


^  +  u  ^  +  u  ^ 
St    1  Sx     2  Sy 


27r, 


=  -  1  - 


ay 


Here  u,  and  Up  correspond  respectively  to  the  horizontal  and  vertical 
components  of  velocity;  while  t  and  tt-.  are  respectively  proportional 
to  the  time  and  pressure.   The  boundary  conditions  which  must  be 
satisfied  at  the  free  surface 

y  =  f-|_(x,T)  , 

where  we  take  the  pressure  to  be  zero,  are 


(2.55) 
and 

(2.56) 


Sf,   Sf-, 

Uo  =  U,  ^ + 

2    1  dx    St 


7r^(x,f-j^,T)  =  0 
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The  boundary  condition  for  the  bottom  is 


(2.37) 


U2(x,0)  =  0 


This  set  of  partial  differential  equations  and  boundary- 
conditions  are  satisfied  by  the  quantities 


^1  "  %^^^  " 


c+  v^(y) 


o 


'gh 


(2.38) 


Ug  =  0  ,    f-^  =  1   , 


[      TT-^     =     1-   Y     , 


where  c  is  a  constant  and  v  (y)  is  a  continuous  non-negative 
function.   These  quantities  define  a  steady  parallel  motion  in  the 
channel  with  the  dimensionless  depth  of  the  liquid  equal  to  unity; 
and  we  will  refer  to  this  flow  as  the  equilibrium  flow.   The 
velocity  v  (y)  gives  the  shear  in  the  axial  velocity  component  and 
it  is  also  a  measure  of  the  departure  of  the  flow  from  a  uniform 
state  defined  by  the  velocity  c. 

Let  us  proceed  to  linearize  the  equations  (2. 33)- (2. 37)  with 
respect  to  the  flow  defined  by  (2.38).   That  is,  let  us  write 


(2.39) 


f  u^  =   u^Cy)  +  u(x,y) 


u^  =  v(x,y) 


TT^  =  1-y  +  M^.-y) 


f^  =  1  +  f{x) 
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and  assume  that  each  of  u,  v,  ^   and  f  Is  Independent  of  the  time 
and  small  with  a  common  order  of  magnitude.   Let  us  substitute 
these  quantities  in  (2. 33)-(2. 37)  and  retain  only  first  order  terms. 
The  result  of  the  linearization  of  equations  (2. 33)-(  2. 34)  is 


(2.40)  ^x'^'^y  =  °  ' 


(2.41)  uu  +u  V  =    -   h      , 

^  '  o  X    oy      ^x  ' 


(2.42)  u^v^        =.  -  ^y 


The  condition  to  be  satisfied  at  the  bottom  of  the  channel  is 

(2.43)  v(x,0)  =  0  . 

Subject  to  the  linearization,  the  free  surface  conditions  become 
conditions  to  be  satisfied  at  y  =  1.   In  place  of  (2.35)  we  have 


(2.44)  v(x,l)  =  u^(l)f^(l,x)  , 


and  from  (2.36)  we  have 

-f  +<|)(x,l  +  f )  =  0 

which  after  differentiation  and  removal  of  higher  order  terms 
becomes 


(2.45)  -  f^(x)  +<f)^(x,l)  =  0 
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The  elimination  of  u,  v  and  f  from  the  linear  equations 
(2.40)-(2.45)  shows  that  (t)(x,y)  must  satisfy 

(2.46)    A   1_  ^  ^(x,y)  +   l—  ^  =  0  ,     -00  <  X  <  CD  , 
5y  u^(y)  ^y  u^(y)  b^ 


0  <  y  <  1 


subject  to  the  boundary  conditions 


(2.47)  (t.y(x,l)  +u^(l)^^^(x,l)  =  0  , 

and 


(2.48)  <t)y(x,0)  =  0 


These  equations,  (2. 46)-(2. 48) ,  constitute  a  particular  case  of  the 
system  (2.1) -(2. 3)  which  we  have  already  analyzed. 

The  system  (2. 46)-( 2. 48)  is  satisfied  by  ^{x,y)   =  const. 
This,  according  to  the  linear  theory,  is  a  necessary  and  sufficient 
condition  for  the  existence  of  an  equilibrium  flow  in  the  channel. 
We  now  ask:   What  conditions  must  u  (y)  satisfy  in  order  to  insure 
that  (j)  =  const,  is  the  only  bounded  solution  of  (2. 46)-(2. 48)? 
As  we  have  seen,  the  answer  to  this  question  can  be  made  to  depend 
on  the  disposition  of  the  eigenvalues  of 


(2.50)  ^y(l)  -  u^(l)A^(l)  , 

(2.51)  ^y(O)  =  0  . 
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It  Is  easy  to  see  that  all  of  the  eigenvalues  of  (2.49)- 
(2.51)  are  real  and  that  A  =  0  is  always  an  eigenvalue  for  which 
we  can  take  ijj     =  1  as  the  corresponding  eigenfunction.   It  is  also 
evident  that  an  eigenfunction  corresponding  to  a  positive  eigen- 
value cannot  vanish  at  the  end  points  y  =  0,  1  nor  can  it  vanish 
at  y  =  Ti  where  0  <  y  =  ri  <  1.   If  it  did,  that  is,  if  ipi-r])    =  0  we 
would  have 


f  *'^'l?:T7-r^*"*^ 


0 


<M 


f  ^^y 


0  %'='' 


which  after  integration  and  use  of  the  boundary  conditions  gives 


\7     dy  =  A  r  yzl  dy  , 

0     ^o^y)  0  %^y^ 


a  contradiction  if  A  is  positive.   (All  of  the  elgenf unctions, 
except  for  a  multiplicative  factor,  must  be  real.)   It  follows  that 
if  A  is  positive,  we  can  divide  (2.49)  by  f   and  integrate  so  as  to 
get 


0 

^0 

0    ^0 

0 

y 
2,2 

dy              =  A  f   ^     , 
^0    ^0 

A 

r       1 

/ 

L  0 

dy 

1 

=   f    4^dy 

0     <^ 
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which  shows  that  a  positive  eigenvalue  cannot  exist  if 


(2.52) 


/i|<i 


0        o 


Furthermore,    if  we   integrate    (2.49)   we   find 


or 


u^(y)     y 


o 


0 


0        o 


0 


(2.53) 


^(1) 


-f 


■^  dy 


0        o 


=   0 


This  can  be  regarded  as  an  equation  for  the  determination  of  the 
eigenvalues  of  (2,  49)-(2.  5I ) .   If  we  expand  ip   into  a  power  series 
in  A  and  substitute  in  (2.53)  we  find  that 


(2.5^) 


0 


dy  _ 
~?  - 

u 
o 


is  the  condition  for  a  multiple  eigenvalue  at  A  =  0.   That  is,  if 


J 


dy 

0     O 


-  1 


is  negative  and  then  made  to  Increase  to  zero  by  changing  u  ,  the 
negative  eigenvalue  of  least  absolute  magnitude  must  approach  zero, 
It  follows  from  the  above  remarks  and  the  criterion  deduced 
from  (2.32)  that  ^   =   const,  is  the  only  bounded  solution  of  (2.46)- 
(2.48)  if 
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u   =  

o 


satisfies 

1 

(2.55)  gh/ ^f-TT^ll  • 

0   t'^+^o^y^^ 


If  V  (y)  =  0  the  condition  (2.55)  shows  that  the  linear 
hydrodynamical  theory  predicts  that  the  uniform  flow  defined  by 


c 

u  =  


/gh 
V  =  0 

is  a  unique  bounded  flow  if  the  speed  c  is  not  less  than  the 
critical  speed  v/gh  where  g  is  the  acceleration  due  to  gravity  and 
h  is  the  depth  of  the  liquid. 

Another  analysis  of  the  foregoing  channel  problem  and  the 
problem  for  a  liquid  with  non-constant  density  can  be  found  in  a 
report  by  Peters  [9]  which  contains  a  more  detailed  discussion  of 
(2-55)  and  other  results. 


3.  Churchill ' s  Method 

If  certain  conditions  are  satisfied,  the  eigenvalue  problem 
{2A)-{2.6) ,    namely 

(5.1)     ^  P(y)  ^  ^(y)  +  q(y)^  -Ar(y)V/  =  0  ,      0  <  y  <  1  , 


2k 


(3.2) 
(5.3) 


^y(O)  +Pq^(0)  =  a^A^(O)  , 


^y(l)  +  P^^(l)  =  a^At(l) 


can  be  reduced  to  an  ordinary  Sturm-Llouville  problem  by  using  an 
appropriate  substitution  for  ip^y).      For  the  case  where  the  a's  and 
P's  are  real  with  ^-,    ^  0   while  a  >•  0,  Churchill  devised  a  pro- 
cedure for  reducing  (3.1)-(3.3)  to  an  eigenvalue  problem  in  which 
the  eigenparameter  does  not  appear  in  the  boundary  conditions.   We 
proceed  to  show  that  with  a  and  a.,  subject  to  no  restrictions  with 
respect  to  sign,  Churchill's  method  can  still  be  used;  provided 
there  exists  a  positive  function  f    (y)  which  satisfies 


(3.4)      ^p(y)  ^^^(y)  +  q^^-A^r^^  =  0  , 


0  <  y  <  1  , 


(3.5) 


i/    (0)  +  P  ^(0)  =   a  J\   Tp    iO)    , 
o  ^o        o  o  o   '  ' 


(3.6) 

and 

(3.7) 


^^(1)  +  P^^(1)  =  a^A^^^(l)  , 


^^{y)  >   0  , 


0  <  y  <  1 


Suppose  that  A   is  real  and  that  it  is  the  least  non- 
negative  value  corresponding  to  which  such  a  ^  (y)  exists.   Define 


w(y)  by 
(3.8) 


^(y)  =  ■^o(y) 


y 


r    w(Tl)dTl 


+  a 


L  0 


^o(t])p(ii) 
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Under  this  substitution,  with  the  eigenvalue  problem  (  5- 1 )- (  5- 3) ;. 
with  A  7^  A  ,  is  changed  Into  the  one  defined  by 


(3-9) 

(5.10) 


d   1   dw  .  (^  _^  )   w  ^ 
dy  ;:^  dy   ^    °  ^ 
^o  o 


r(0)w(0)  =  a^p(0)w'(0)  , 


(3.11) 


r(l)w(l)  =  a^p(l)w'(l) 


Conversely,  this  standard  Sturm-Llouvllle  problem  Is  transformed 
by  (3.8)  into  the  problem  of  finding  f   such  that 


(3.12)     ^  P^'  +q^- Ar^  =  r^^ 


w'(0)_  _  (,_xja 


r(0)T//^(0) 


(3.13)      ^'(0)  +  6t(0)  =  a  A^(0)+  ^^(0) 


and 


o 


o 


T^ (A  -  ^^]a 


r(0)^^(0) 


(3.1^)     t'(l)+P^t(l)  =  a^A^(l)  +  to(l) 


^'^|L  _  (A-Aja 


r(0)^Q(0) 


If  A  =  A   is  not  an  eigenvalue  of  (5.9)-(3.11)  the  constant  a  can 
be  fixed  by 


(3.15) 


w'(0) 


r(0)^Q(0) 


(A  -  A  )a  =  0 


and  then  (3.9)-(3.11)  is  equivalent  to  (3.1)-(3.3) 
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Suppose  next  that  A  =  A  is  an  eigenvalue  of  (5. 9) -(5-11) 
and  that  the  corresponding  elgenfunctlon  Is  w  (y)  so  that 


(3.16) 


dy  r^2  dy 
^o 


=  0  , 


(5.17) 


r{0)vj^{0)   =  a^p{0)w^(0)  , 


{3.18) 


r(l)w^(l)  =  a^p{l)w^(l) 


The  function 


(3.19) 


r  y 


w^(y)  -  k 


/   r(Ti)^^(Ti)dT]  +a^p(0)^^(0) 


-  0 


w'(0) 

o^  ' 


r(0)^^(0) 


satisfies  (3.16)  and  (5. 17).   The  condition  {3.18)  then  shows  that 
A  is  an  eigenvalue  if  and  only  if 

1 
(3.20)     f    r(y)^^{y)dy-a^p{l)t^{l)  +  a^p(0)^^(0)  =  0 


or 


(3.21) 


Q[i^^{Y),fjY)]    =   0 


if  we  use  the  symbol  defined  in  Section  2.   The  function  T^(y)  which 
corresponds  to  w  is,  in  accordance  with  (5'8)>  defined  by 
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dy  ^^^7^  "  p(y)^^(y)  ' 


and  it  is  easy  to  verify  that  t}/   must  satisfy 

(3.22)  ^  p(y)  ^  ^'(y)  +  q?  -A^r?^  =  kr^^  , 

(3.23)  r  (0)  +  pQ?(0)  =  a^A^5^(0)  +a^k^^(0)  , 

where  k  7^  0.   The  function  ^(y)  is  not  an  eigenfunction  of  (3.I)- 
(3.3).   However,  if  the  set  of  eigenfunctlons  of  (3-l)-(3-3)  is 
extended  by  the  addition  of  '^(y)  then,  as  we  shall  see,  an 
arbitrary  function  can  be  expanded  with  respect  to  the  enlarged 
set. 

Under  our  assumptions,  the  eigenvalues  of  (3.9)-{3-ll)  are 
all  real.   These  eigenvalues  can  be  ordered  with  respect  to  absolute 
magnitude,  that  is, 

and  it  is  well  known  that  the  corresponding  set  of  eigenfunctlons 
is  complete.   If 


(3.25)  Q[^o(y),^o(y)l  ^  0  , 


the  set  of  eigenfunctlons  can  be  designated  by 
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{w^(y)} 


n  =    1,2,5, 


The  eigenfunctions  satisfy 


^  w^(y)w^(y) 


m 


0  p(y)to(y) 


dy  -  0  , 


m  7^  n 


If  F(y)  Is  dlf ferentlable  we  have 


r 


F(t)w^(t) 


dt 


F(y)  = 


00   0  P(t)^o(t) 
ri=l 


f     <(*' 


w^(y) 


0  p(t)tj(t) 


dt 


Hence  If  f(y)  Is  an  arbitrary  dlf ferentlable  function  we  have 


(3.26) 


I       dt 


0 


f(t) 


y 


irm 


dt 


^f 


0  p{t)^;(t) 


\y~  7  w  (t)>  dt 


where 


(3.27) 


7   = 
'n 


/ 

0 


w„(t)  ^  (f  )dt 

o 


0 


^  w^(t) 
n^  ^ 


n  =  1,2, 


dt 


P^, 


o 


In  terms  of  the  set,  {^^(y)},  of  eigenfunctions  of  (5-l)-(3-3), 


n 


Including  f^[Y) ,    the  expansion  (3.26)  Implies 


(3.28) 


OD 


f(y)  =ZZ  ^n^n^y) 

n=0 


where 
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(3.29) 


If 


{3.30)  Q[V'o(t),^^(t)] 


=  f 


r(t)V/^(t)dt-a^p(l)^^(l)+a^p(0)^^(0) 


0 


=  0 


the  elgenfunctlons  of  (3-9)-(5-ll)  can  be  designated  by 


(w^Cy)} 


n  =  0,1,2, . . . 


For  this  case  we  have 


{3.31) 


y 


d 

dT 


0 


f(t) 


y 


j;m 


00 


dt 


f    K flZ7w(t)?dt 

P(t)<(t)  tfery    ^  ^      J 


0 


where 


(3.32) 


7 


r 

0 


w    (t)   A   (-^)dt 
n^    -*    dt    ^^  ^ 
'^o 


n 


J~~^ 


n  =   0,1,2, 


0     P(t)^Q(t) 


dt 


This  can  be  expressed  In  terms  of  the  elgenfunctlons,  [^  (y)),  of 
(3'l)-(3-3)  and  the  generalized  elgenfunctlon  ip[j) .      A  calculation 
shows  that 


(3.33) 


00 


f(y)  =  s^(y)+  s^^(y)  +  5^  s^n^y^ 

n=l 


where 
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Q[f,^n^ 

n 

Q[^o,^] 

(3.3^) 


(3.35) 

and 

(3.36)         s  =  2 ^  ^  ° . 

The  expansions  (5.28)  and  (3-33)  are  particular  cases  of  the 
basic  general  expansion  (2.19).   It  appears  that  the  assumption 
with  respect  to  the  existence  of  ^  (y)  implies  that,  excepting  a 
possible  multiple  pole  at  z  =  A  ,  G(t,y,z)  possesses  only  simple 


poles. 


When  ^  (y)  exists  the  question  of  the  existence  of  a  non- 


trivial  bounded  solution  of  (2.1)- (2. 3)  reduces  to  finding  the 
disposition  of  the  eigenvalues  of  (3.  9)- (3- H)  •   Under  certain 
circumstances  this  disposition  can  be  exposed  more  clearly  by 
integrating  (3.9).   When  A  ^  A   integration  of  (3-9)  namely, 

^       1   "^w    r-x   >  ^   w    d  f    1     dw    ,,   -V  >  ^  ] 


gives 


r^Q  ^  P^o    -^  <-  rf^ 


7^  dy       o^  ?r 
o 


or 


o 
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Suppose  now  that  if  A  Is  a  positive  eigenvalue  such  that  A  >  A 


o 


then  the  corresponding  eigenfunction  does  not  vanish  in  the  interval 
0  j^  y  _^  1.   When  this  condition  prevails  we  can  multiply  (3.37)  by 
rf  /f   to  obtain 


o' 


— i-  -3—  p^     -3—    (-r-)    =    (A  -  A    )rf 
Tp    dy  ^'^o   dy    ^1^    '         ^  o"^    ^o 


Another  integration  here,    produces 


(3.38) 


0 


^^  2r 


(/) 


_d_   ,  t  , 
_dy    ^^^_ 


-i2 


P^^dy 


=    (^-^o) 


][    r^^dy-  a^p(l)i//^(l)  +  a^p{0)^^(0) 


This  shows  that,  subject  to  the  supposition  above,  there  cannot  be 

a  real  eigenvalue  A  such  that  A  >  A  if 
°  o 


Q[^.,^^]  =  /   rt^dy-  a  p{l)^^(l)+'a^p(0)^?{0) 


o'  '  o 


o 


o'- 


0 


is  negative.   When  this  criterion  is  applied  to  (2. 49)-(2. 5I)  it 
yields  (2.52). 


32 


4.  Application  of  a  Transform  Method 

The  system  which  was  analyzed  In  Section  2,  namely 


(4.1)   |-  p(y)  A  (|)(x,y)  +  q(y)^  +  r(y)  ^  =  0  , 
By      ay  ^^^ 


-OD   *=  X  <  00   , 

0  <  y  <  1   , 


(4.2)  <t)y(x,0)  +  a^(t)^(x,0)+  P^(t)(x,0)  =  0  , 

(4.3)  (t)y(x,l)  +  a^(t)^^(x,l)  +  P-^())(x,l)  =  0  ; 

can  also  be  studied  by  using  a  generalized  Fourier  transform  method. 
Let  the  right-hand  Fourier  transform  of  <f)(x,y)  be 

00 


$(>^,y)  =  f    e^^^<f)(x,y)dx 


0 

where  Im  A  =  a  >  0.   Let  the  left-hand  transform  of  <j)  be 

0 


f   (A,y)  =   f    e^^^Hx,y)dx 


-00 


where  Im  A  =  b  *=  0.   By  taking  the  magnitudes  of  a  and  b 
sufficiently  large  these  transforms  exist  for  a  twice  dlf ferentlable 
function  ({)  of  any  exponential  order.   The  recovery  formula  for 
(j){x,y)  Is 

oo  +la  00  +lb 

-00  +la  -oo  +lb 
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If  ({)  and  its  first  derivative  with  respect  to  x  are  of  exponential 
order,  the  application  of  the  right-hand  transform  to  (4.1)- (4. 3) 
for  X  >  0  shows  that  ^  must  satisfy 


(^•4)   4-  P(y)ct)  O^.y)  +q?-A  r^  =   r[(t)  (0,y)-  iA(f)(0,y)] 
C3y      — y         ~     —       ->^ 


for  0  ^  y  "*  1  subject  to  the  boundary  conditions 

(4.5)  $^(A,0)  +p^(A,0)  =   CL^-^'^{^,0)  +a^[^JO,0)  -  imo,0)] 
and 

(4.6)  ^  (A,1)  +  P^J(A,1)  =  a^A^(A,l)  +  a^[(|)^(0,l)-  lA(f)(0,l)]  . 

Similarly,  the  application  of  the  left-hand  transform  to  (4.1)-(4.3) 
for  X   <   0   shows  that  J  must  satisfy 


(^.7)   ^ 


-1 
2 


^  P$   +qJ^-A  rf^  =  -r[(t)^(0,y)  -lA(t)(0,y)] 

for  0  "^  y  <  1  with  the  boundary  conditions 

(4.8)  ^   (A,0)  +P^  (A,0)  =  a  A^^  (A,0)  -  a  [(|)  (0,0)  -iA(t)(0,0)] 
~ly  '-*"^1  °   ~1  OX 

(4.9)  $   (A,l)  +p^^(A,0)  =  a^A^  (A,l)  -a^[(l)^(0,l)  -  1A(|)(0,1)] 


From  the  equations  and  boundary  conditions  which  "^  and  '^ 
must  satisfy  it  is  evident  that 

I^(A,y)  =  $(A,y)  . 
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Therefore  we  see  that 


(^•10) 


H^,y)    -    -  ^   Te-^^^  ©A,y)dA 


M 


where  M  is  the  path  M  =  M-,  +  Mp  shown  In  Fig.  4.1.   The  lines  M-,  and 


M. 


<y 


X-plane 


M. 


Fig.  4.1 

Mp  are  parallel  to  the  real  axis  of  the  A-plane  and  their  respective 
distances  a  and  |b|  from  the  real  axis  can  be  adjusted  to  admit 
functions  cj)  of  various  exponential  orders.   Also,  by  proper  choice 
of  a  and  b  we  can  allow  the  exponential  order  of  behavior  of  ^   as 
X  — >  00  to  be  different  from  that  as  x  — >  -oo  . 

It  can  be  shown  that  ^  is  expressible  as  a  ratio 


?(A,y) 


^(y,A) 

a)(A) 


in  which  each  of  Tp   and  ca   is  an  entire  function  of  the  complex 
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variable  A.  The  zeros  of  a)(A)  are  just  the  eigenvalues  of 

(4.11)  ^  P^y(y,A)  +q^-  A^rt  =  0  ,  0<y<l, 

(4.12)  V'y(0,A)  +p^^(0,A)  =  a^A^^(0,A)  , 

(4.13)  ^y(l,A)  +P-^^(1,A)  =  a-|^A^^(l,A)  . 


If  the  disposition  of  these  eigenvalues  is  known  and  if  the  be- 
havior of  ■'/'(y, A)/(d( A)  as  A  — >  00  can  be  estimated,  then  the  be- 
havior of  ^{Xf-y)    with  respect  to  x  can  be  found  from  (4.10)  by  using 
the  theory  of  residues.   Hence  our  problem  is  again  reduced  to  a 
study  of  the  elgenfunctlon  system  which  was  introduced  in  Section  2. 
Suppose  for  example  that  the  parameters  of  (4. 11)- (4. 13)  are 
such  that  all  of  the  eigenvalues  are  pure  Imaginaries  not  including 
A  =  0.   For  a  bounded  solution  ^   we  can  take  a  =  |b|  =  e  so  small 
that  M  does  not  contain  any  of  these  eigenvalues.   Then,  since  our 
assumptions  on  the  order  of  (j)  and  its  first  derivative  with  respect 
to  X  imply 

e 
lim    f     e-l(N+iTi)x  ^(N+lTi,y)dTl  =  0  , 
N — >±oo 


-e 


it  follows  that 


H^>y)    =    -  ^   1^    e"^^^  $(A,y)dA  =  0 
M 


is  the  only  bounded  solution  of  (4.1)-(4.3)  if  (4 . 11 )-(4 . 13) 
possesses  no  real  eigenvalues. 
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It  should  be  noted  that  the  above  transform  method  does 
not  require  a  knowledge  of  the  completeness  or  lack  of  completeness 
of  the  elgenf unctions  of  ( 4, ll)-( 4, I3) .   Hence  the  infinite  trans- 
form method  appears  to  be  a  direct  one  for  at  least  the  investiga- 
tion of  bounded  solutions  of  (4.1) -(4. 3).   However,  this  method 
starts  with  assumptions  about  the  behavior  of  ^   and  its  derivatives 
at  infinity  whereas  the  generalized  eigenfunction  methods  of 
Sections  2  and  5  do  not  require  such  initial  assumptions.   Further- 
more, the  transform  method  depends  on  an  estimate  of  the  transform 
^(A,y)  as  A  — >  (x>  .      A  consideration  of  these  facts  and  a  comparison 
of  a  transform  method  with  the  method  of  Section  2  as  these  methods 
are  conceptually  applied  to  problems  involving  domains  more  general 
than  the  doubly  infinite  strip  domain,  leads  to  the  conclusion  that 
the  eigenfunction  method  is  more  fundamental.   For  example,  the 
eigenfunction  method  can  be  applied  to  problems  involving  a  rec- 
tangular domain  for  which  a  finite  transform  method  is  not  directly 
effective  because  an  inversion  formula  for  the  transform  is  not  at 
hand  and  needs  to  be  derived.   The  derivation  of  such  a  formula  is 
in  general  equivalent  to  establishing  an  eigenfunction  expansion. 
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